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We study the renormalization group flow of the interactions in the two-dimensional t-t' Hubbard 
model near half filling in a iV-patch representation of the whole Fermi surface. Starting from 
weak to intermediate couplings the flows are to strong coupling with different character depending 
on the choice of parameters. In a large parameter region elastic Umklapp scatterings drive an 
instability which on parts of the Fermi surface exhibits the key signatures of an insulating spin 
liquid (ISL), as proposed by Furukawa et al., rather than a conventional symmetry-broken state. 
The ISL is characterized by both strong d-wave pairing and antiferromagnetic correlations, however 
it is insulating due to the vanishing local charge compressibility and a spin liquid because of the 
spin gap arising from the pairing correlations. We find that the ISL is a consequence of a Fermi 
surface close to the saddle points at the Brillouin zone boundaries which provides an intrinsic and 
mutually reinforcing coupling between pairing and Umklapp channels. 



I. INTRODUCTION 

The Landau theory is widely used to describe Fermi 
liquids even when the interactions are strong, but it 
cannot be justified 'a priori'. The cuprate high-T c su- 
perconductors show clear deviations from Landau the- 
ory in the normal state and it has long been argued 
that the key to understanding these materials lies in 
the breakdown of Landau theorycl. One possible cause 
is a symmetry-breaking instability such as magnetic or- 
der. But in experiments on underdoped cupratescl, the 
marked deviations from Landau theory, such as the on- 
set of the spin gap and gaps in the ARPES spectra 
near the saddle points of the Fermi surface (FS), appear 
without an obvious symmetry-breaking. This raises the 
question whether a breakdown of Landau theory with- 
out symmetry-breaking is possible. Actually one exam- 
ple is well known and understood, the insulating spin 
liquid states of even-leg ladder systems at half-filling, 
which have only short range magnetic order and unbro- 
ken translational symmetrytjcl. The key to this behavior 
are elastic Umklapp scattering processes across the FS 
which open up a charge gap at half-filling in addition to 
the spin gap caused by the pairing instability. In this 
paper the role of these processes in a two-dimensional 
system will be carefully examined. 

Renormalization group (RG) methods allow an analyt- 
ical treatment and, although the one-loop approximation 
is in principle applicable only at weak coupling, we can 
hope to learn about possible instabilities at the strong 
to intermediate couplings that apply in the cuprates. 
Such methods have long been successfulljy.-applied to one- 
dimensional models. The first attemptsLTEl to extend this 
analysis to two dimensions were made shortly after the 
discovery of high temperature superconductivity. They 
focussed on the dominant role of scattering processes in- 



volving the Fermi surface regions in the vicinity of van 
Hove singularities. 

Limiting the two-dimensional FS to just two patches 
reduces the problem to the flow of a small number of cou- 
pling constants which can be handled analytically. For 
repulsive interactions there are two possible fixed points 
involving flows either to weak coupling or to strong cou- 
pling. The possible relevance of the latter to the,«iprates 
was emphasized by three of the present authors^. They 
showed that under certain conditions the local charge 
compressibility flowed towards zero, indicating that here 
too Umklapp scattering opened up a local charge gap. 

A proper treatment requires that the flow of interac- 
tions involving the whole two-dimensional Fermi surface 
be included. Already several RG investigations using a 
discretization of the Fermi surface into N patches with 
N < 32 have been made. Zanchi and SchulzEil studied 
the RG-fiows of a 32-patch weak coupling Hubbard model 
with only nearest-neighbor (n.n.) hopping in the kinetic 
energy term. They found a crossover between an antifer- 
romagnetic (AF) ordered groundstate to a d^a^s -paired 
superconducting (SC) groundstate as the electron density 
was lowered away from half-filling. .Recent more exten- 
sive results by Halboth and Metznenl3 have largely con- 
firmed the Zanchi-Schulz results, extending them to the 
case where there is a small next-nearest-neighbor (n.n.n.) 
hopping as well and investigating possible incommensu- 
rate AF-orderings. Although in both these investigations 
Umklapp scattering was included, the possibility of a 
fixed point behavior which would be similar to that of 
the two-leg ladder was not explicitly considered. 

In this paper we will use a one-loop RG method with a 
discretization of the FS into A-patches {N = 32 — 96) to 
examine the flow of the coupling constants and suscep- 
tibilities under various starting conditions. Throughout 
we take a substantial value for the n.n.n. hopping ampli- 
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tude, t'. On the one hand this is a realistic value for the 
cuprates. Secondly it moves the critical density, where 
the saddle points are at the FS, away from half- filling 
so that the saddle point effects are not mixed with nest- 
ing effects on the zone diagonals, as occurs when one 
sets t' = 0. When t' is substantial, one can distinguish 
three density regions. The simplest is the strongly doped 
region where the saddle points lie above the Fermi en- 
ergy and Umklapp scattering is unimportant. Here the 
leading instability is to d-wave SC — a form of Kohn- 
Luttinger instability, in agreement with previous studies. 
We call this the d-wave dominated regime. A second rela- 
tively straightforward density regime is the weak doping 
regime close to half-filling, where the approximate nest- 
ing of FS segments near the zone diagonals dominates 
and a AF instability is favored — again in agreement 
with previous studies. We call this region the approxi- 
mate nesting regime. 

The intermediate regime is most interesting and will be 
the focus of this work. In this case the saddle points lie 
slightly below the Fermi energy and Umklapp processes 
involving these FS regions are highly relevant. We call 
this density region the saddle point regime. As in the 
case of the two-leg ladder these Umklapp processes act 
to reinforce d-wave pairing so that this channel competes 
strongly with the AF instability. If one looks only at 
these two instabilities in the one-loop RG, it is not pos- 
sible to decide which dominates. In the case of the two- 
leg ladder the uniform (Pauli) spin susceptibility flows 
to zero indicating that pairing instability prevails — a 
result confirmed when bosonization methods are used to 
examine the strong coupling state below the critical scale 
in the one-loop scheme. Similarly we find in the present 
case that an examination of the uniform susceptibility fa- 
vors an assignment of the strong coupling fixed point to 
the class of the two-leg ladder. Further the local charge 
compressibility defined for these FS segments appears to 
scale also to zero. Below we present a detailed examina- 
tion of this saddle point regime. We argue that rather 
than the simple crossover between d-wave SC and AF or- 
der as the density varies, found by previous authors, an 
interpretation in terms of the formation of an insulating 
spin liquid (ISL) which truncates the FS segments near 
the saddle points is justified. This ISL can be viewed as a 
form of d-wave RVB (Resonance Valence Bond) state as 
in the case of the two-leg ladder. Such a state represents a 
clear violation of the Landau theory, which does not rely 
upon a translational symmetry breaking mechanism. 

Clearly an instability that partially truncates the 
Fermi surface with a charge gap can be seen as a fore- 
runner of the Mott insulating state which occurs for in- 
termediate to strong interactions. Since our motivation 
is to understand better the phase diagram of the high- 
T c cuprates, we are most interested in such instabilities. 
However we are aware that there are other instabili- 
ties which appear in a weak coupling theory driven by 
the diverging density of states (DOS) at the van Hove 
points. These are the Stoner instability to ferromag- 



netism and, as remarked by Halboth and Metznerliia, the 
Labbe-Friedel or Pomeranchuck instability from square 
to rectangular symmetry. These split the saddle points 
when the Fermi energy lies near the van Hove singular- 
ity. This is a drawback of using a weak coupling approach 
to describe an intermediate to strong coupling problem. 
We will simply ignore these DOS-related instabilities in 
the forward scattering channel and concentrate on those 
which we believe are more relevant as weak coupling sig- 
natures of the intermediate to strong coupling problem. 

Finally, a defect of our one-loop approximation is 
that it does not lead to a description of the strong cou- 
pling phase of the system which it predicts. There are 
many approaches in the literature which attempt to con- 
struct a theory of such a state that|-we_can loosely call 
a lightly doped d-wave RVB state£lla't3~t3. Our aim here 
is rather different and seeks to complement these strong 
coupling theories by examining the approach from the 
strongly overdoped regime which behaves as a conven- 
tional Landau-Fermi liquid with an instability towards 
weak-coupling d-wave superconductivity. The question 
we seek to address is the form of the instability in a 
Landau-Fermi liquid which leads to this doped RVB 
state. How does it differ from a simple d-wave super- 
conducting instability and how does the proximity to the 
Mott insulating state at half-filling manifest itself ? 

II. THE MODEL AND ITS FERMI SURFACE 

The kinetic energy of the t — t' Hubbard model is given 
by the tight-binding dispersion 

e(k) = — 2£(cos k x + cos k y ) + At' cos k x cos k y — /i (1) 

with nearest neighbor (n.n.) hopping t, next nearest 
neighbor (n.n.n.) hopping t' and chemical potential /i. 
Typically we choose t' — 0.3t which yields a small convex 
curvature of the FS around (tt,tt) at higher fillings. An- 
other essential curve is the Umklapp surface (US) which 
connects the van Hove points with straight lines. If the 
FS crosses this line, two particles at the FS can be scat- 
tered from one side of the US to the opposite one in an 
elastic process. As we will see, these additional scatter- 
ing channels then enhance the scale of the transition to 
a strong coupling regime. 

The initial interaction is taken to be a simple on-site 
repulsion 

H u = u z2 n S,t n x,l' ( 2 ) 

X 

which is constant in /c-space. The effective interaction 
will develop a pronounced /c-space structure in the RG 
flow. 

In recent years there have been several RG approaches, 
to the 2D Hubbard model. Schulztl and Lederer et al.H 
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studied the RG flow of the processes connecting the sad- 
dle points emphasizing the divergence of both AF and d- 
wave pairing correlations. DzyaloshinkiiQ discussed the 
weak coupling non- Fermi liquid fixed point of such .-a 
model. Similar studies have been given by Alvarez et al.t£l 
and Gonzalez et al.lHI. Later on, in a relatad-,formalism 
based on parquet equations, Zheleznyak et alJ13 examined 
the interplay between critical scales and effects of the FS 
curvature for a quasi 2D model restricted to approxi- 
mately flat FS faces close to half-filling. Another study 
of nesting effects between flat-FS segments was gijron by 
Vistulo de Abreu and Doucott3. Zanchi and SchulztU pre- 
sented the first fully two-dimensional treatment, based 
on Polchinski's RG equation. They studied the 2D Hub- 
bard model with t' = and found two different regimes 
with dominant AF in the one and d-wave pairing correla- 
tions in the other. A more detailed analysis of th&leading 
instabilities was given by Halboth and Metznerlia using 
RG equations for Wick ordered functions. In this pa- 
per we study the RG flow for the one-particle irreducible 
(1PI) vertex functions for the model given by (|j) and 
(||), to investigate the possibility of a novel strong cou- 
pling phase which is a precursor of the Mott insulating 
state^-as suggested by the two-patch study by Furukawa 
et al.E2l. A self-contained and detailed discussion of the 
RG technique we use here is given elsewherec£l. 

III. THE TWO-PATCH MODEL REVISITED 

We start with a brief discussion of the dominant mecha- 
nisms for the case where the FS is at the saddle ppin^. 
These are most transparent in the two-patch mode£3 : B : EI, 
where only small phase space patches around the saddle 
points at (n, 0) and (0,7r) are kept. Neglecting possible 
frequency dependence we can approximately describe the 
scattering processes within and between the two patches 
by four coupling constants, gi ■ ■ ■ 34, depicted in Fig. [l| 
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FIG. 1. The relevant scattering processes in the two-patch 

model. The gray semi-circles denote the phase space patches 
around the saddle points. The interactions are assumed to 
be spin-independent and constant over the patches. In this 
notation the spin of the initial and final particle connected by 
an arrow has to be the same. 

The main terms which drive the one-loop RG flow 
of these vertices are: a) the particle-particle loop do 
with zero total incoming momentum, which diverges like 
log 2 (A /A) with decreasing energy scale A < A due to 



the van Hove singularity in the density of states and b) 
the particle-hole loop with momentum transfer (-7T, it) de- 
noted by d\ which, in presence of a small but nonzero t' , 
diverges like log(Ao/A) with a large prefactoiO. Keep- 
ing only these two contributions, and denoting y — 
log(Ao/A), so that decreasing A means increasing y, we 
obtain the RG flow equations 

.gi = 2d 1 g 1 (g 2 - gi) , (3) 

92 = d\ (gi + gi) , (4) 

53 = -2d .9334 + 2d'i g 3 (2g 2 - g%) , (5) 

.94 = -do [gj + g\) . (6) 

where gi = dgi/dy, and do, d\ > 0. It is useful to briefly 
review_.the analysis of the two-patch model by Furukawa 
et al.E3. The second term on the right hand side of (|^) 
enhances the basin of attraction of the strong coupling 
fixed point. Starting from the onsite repulsion g\—gi — 
.93 = gi~U given by (§) , the coupling constants diverge 
at a scale A c : g% — > +00, g± — > —00 and 172 — * +00; 
31 diverges more slowly. Initially there is a competition 
between the two terms on the r.h.s. of (||), but the r.h.s. of 
(||) is always negative and thus decreases 54. Eventually, 
54 becomes negative; then both terms in have the 
same sign, which accelerates the flow to strong coupling. 

For incoming and outgoing particles directly at the 
saddle points (^-processes correspond to both Cooper 
and Umklapp processes. However, away from the saddle 
points we can distinguish between Cooper processes with 
approximately zero total incoming momentum driven 
through the particle-particle channel and Umklapp pro- 
cesses with momentum transfer ks (tt, it) driven by the 
corresponding particle-hole channel. ^From this point 
of view, Eq. (g) states that for incoming and outgoing 
wavevectors near the saddle points the Umklapp and 
the d x 2_ y 2-wave Cooper channel are coupled and mutu- 
ally reinforce each other through the 53 and (74-processes 
which belong to both channels, thereby increasing the 
critical scale A c . In fact the divergence of the Umklapp 
scatterings processes implies a divergence of the d-wave 
couplings and vice versa. 

An analysis of the susceptibilities shows a competition 
between divergences in the d x 2_ y 2 -pairing and AF chan- 
nel controlled by the flow of the combinations (173 — 54) 
and (<?2 + 33) respectively. In this case of competing sin- 
gularities it is not clearcut which of them dominates. Fu- 
rukawa et al. proposed to resolve the issue by examining 
the uniform spin susceptibility and the charge compress- 
ibility. For a not-too-weak value of U/t and t'/t they 
found both were driven to zero by the pairing and Umk- 
lapp processes respectively. On this basis they assigned 
the fixed point to be in the same class as that of the re- 
pulsive two-leg ladder at half-filling. In that system, the 
one-loop RG also exhibits competing and equally strong 
divergences in the d-wave pairing and AF channels, but 
the groundstate is known to be an insulating spin liquid 
(ISL) from a bosonization treatment of the strong cou- 
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pling regimeotil. This form of groundstate is already sig- 
naled in the RG calculation by the suppression towards 
zero in the uniform spin susceptibility and the charge 
compressibility. The ISL in the two-leg ladder at half- 
filling is a form of RVB (resonance valence bond) state 
with an approximate d-wave pairing symmetry, but with- 
out any explicit translational or gauge symmetry break- 
ing. 



IV. THE TECHNIQUE 

For the iV-patch analysis we use a Wilson RG flow for 
one-particle irreducible (1PI) vertex functions. The full 
RG flow associates to every energy scale A below the 
bandwidth Ao an effective interaction for the particles 
with energies e(fc) below A, in a way that the generating 
functional for the Green functions remains independent 
of the scale A. Because of this exact invariance, the effec- 
tive interaction is no longer just quartic but an infinite 
power series in the fields. The full RG can be expressed 
as an infinite hierarchy of differential equations for the 
1PI m-point vertex functions. Here we study a trun- 
cation of this infinite system in which only the 2- and 
4-point functions are kept. For a derivation of the full 
flow equations and a discussion of this truncation seeE3. 
Here we just state the results, which are rather simple. 
Because our model is two-dimensional, continuous sym- 
metries cannot be broken by long-range order at any pos- 
itive temperature. Therefore the effective action must 
be gauge-invariant and invariant under spin rotations, 
hence the four-point function is determined by the func- 
tion V\(u>i, ki, u)2, &2) fe) which describes the scatter- 
ing of two incoming particles (wi, fci, oi) and (0*2, &2, °2) 
into two outgoing particles (0^3,^3,(73) and (074, £4, 04) 
where <j\ — 03 and 02 — cr^, lo a — + ^2 — ^3, and k± 
is given by momentum conservation as ki = k\ + ki — k-& 
modulo reciprocal lattice vectors. Because the spin of 
particle 1 (first incoming) is the same as that of particle 
3 (first outgoing), and similarly for 2 and 4, we may draw 
the vertex corresponding to Va as in Fig. |2[ where the 
solid fermion lines going through at the top and the bot- 
tom of the vertex indicate that spin is conserved along 
these lines. 



G A (k,iu n ) 



XA(fc) 



iu> n - e(k) - x A (fc)S A (fc, iu n ) 



(7) 



where XA(fc) = l-{exp[(|e| - A)/(0.05A)] + 1} 1 cuts off 
energies belaw A. The other line stands for a single-scale 
propagatorEj 



S A (k,ico r , 



Xh{k) (iuj n - e(k)) 



(iuj n - e(k) - XA(fc)S A (fc, iuj n )) 2 



(8) 



Because it contains the derivative of the cutoff function 
with respect to A, Sa is nonzero only at energies close 
to A. Since there are two possibilities for assigning Ga 
and Sa to the internal lines, each graph stands for two 
contributions. Apart from that the usual diagrammatic 
rules hold: the graph with the fermion loop gets a factor 
2 from the spin trace and a minus sign. 

The contributions to the self-energy have the graphical 
representation shown in Fig. ^. Here the internal line 
stands for a single-scale propagator S\. 

In the main part of this paper we will neglect selfen- 
ergy corrections to the propagator. Then G(itu n ,k) = 
Xh{k)(ioJ n — and the single scale propagator is 

simply Sa = <9Ga/<5A. In Appendix A we show some 
results for a flow with the real part of the selfenergy on 
the FS taken into account. A more complete study in- 
cluding selfenergy effects, in particular the effects of the 
wave function renormalization, is underway. 




(a) 
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FIG. 3. The contributions to the right-hand side of the 
RGDE. (a) the particle-particle term (b) the crossed parti- 
cle-hole term (c) the direct particle— hole terms; the first of 
these three graphs gets a factor —2 because of the fermion 
loop. 



FIG. 2. The vertex corresponding to Va(/ci, ^2,^3) 

The contributions to the right hand side of Va = ^Va 
can then be represented graphically as in Fig. ||. In these 
graphs, one of the internal lines represents a full electron 
propagator 
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FIG. 4. The contributions to the selfenergy 
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We want to emphasize that our RG method does not 
rely on any form of scale invariance or scaling ansatzes. 
The RG we set up in Ref . ^ provides an exact rewriting 
of the generating functional in terms of the effective ac- 
tion. The approximations we make in the present paper 
are (i) we discard the 1PI m-point functions with m > 6, 
(ii) we project the four-point function to the Fermi sur- 
face and frequency zero (see the next section), and (iii) 
we neglect the selfenergy corrections. We discuss the jus- 
tification of (i) and (ii) in detail ir£3. The justification for 
(ii) is a standard RG argument, and the full momentum 
dependence can in principle be reconstructed by calcu- 
lating susceptibilities and related quantities. The justifi- 
cation of (i) is less trivial but possible for curved Fermi 
surfaces and in a specific scale rangecJ. (iii) is an approx- 
imation on which we shall improve in a further paper. 

In systems with a Cooper instability, the flow always 
tends towards strong coupling at a sufficiently low scale. 
This happens even in repulsive systems because of the 
Kohn-Luttinger effect. However, in repulsive, and ini- 
tially weakly coupled systems, the flow stays in the weak 
coupling regime down to a very low scale which may never 
be reached because the temperature, which acts as a nat- 
ural infrared cutoff, stops the flow before that (in the 
usual Kohn-Luttinger effect, this scale is at most of or- 
der e - const / u ■ se e Refs. p2|]23| ). In this case the system 
stays weakly coupled above a certain temperature, and 
components of the Fermi surface limits of the four-point 
function can be. identified with the Landau interaction 
function f(k, k'f^B. 

When the four-point function flows to strong coupling, 
the critical scale A c where the coupling constants diverge 
gives an estimate for the scale where the quasiparticlcs 
will be strongly modified or entirely destroyed (e.g. for 
the superconducting transition, a gap opens up). 

The general picture we find in this model is that the 
flow always tends to strong coupling, but that for fillings 
where Umklapp scattering is favored by the geometry of 
the Fermi surface, the critical scale is strongly enhanced. 

We note that in our RG method the temperature T is 
retained as a physical parameter and that the RG pro- 
cedure of decreasing the scale A is a priori not related 
to changing the temperature. The four-point vertex at 
scale A is the effective interaction for the modes with 
energy below A, and at the same time it is the four- 
point function with infrared cutoff A. For A <C T, one 
should use the second interpretation. Note that because 
we do not impose a scale-dependent cut-off on the fre- 
quencies the flow does not stop exactly at the point when 
A decreases below the smallest fermionic Matsubara fre- 
quency, u>i = ttT. 

V. NUMERICAL IMPLEMENTATION 

Next we describe the practical implementation of this 
RG scheme for the 2D Hubbard model. First we define a 



phase space discretization following Zanchi and Schulzliil. 
The idea is to discretize the Brillouin zone (BZ) into N 
segments centered around N lines. Each line with index 
k 6 {1...JV} starts from the origin(r-point) in a cer- 
tain angular direction and from the Y point (±tt, ±7r) so 
that the lines meet at the Umklapp surface. All phase 
space space integrations with measure d 2 k/ (2it) 2 are per- 
formed approximately as sums over the lines and integra- 
tions over the radial direction. These imply Jacobians for 
polar coordinates with respect to the T- or the K-point, 
respectively. 




FIG. 5. The Brillouin zone, Fermi and Umklapp surface 
and the lines in the patch centers for N = 32. 

The interaction vertex V\(ki, fe, fe) depends on two in- 
coming wavevectors k\ and k 2 and one outgoing wavevec- 
tor fc 3 lying in segments k\, k 2 and fc 3 respectively (here 
we have already projected the frequencies to zero). The 
fourth wave vector k^ is fixed by momentum conserva- 
tion. In a next approximation we select a large but finite 
number of coupling constants representative for^certain 
regions in the space spanned by k\, k 2 and k^c3. We 
choose to take these wave vectors as the crossing points 
of the lines fci, k 2 or k% with the Fermi surface (FS), 
i.e. k F (ki), k F (k 2 ) and k F (kj,), which lie at the centers 
of the corresponding FS patches. By Taylor expansion 
and power counting argumcntsc2l, the leading part of the 
flow is given by the coupling functions on the Ferroi,sur- 
face and at zero frequency. Thus we approximateE3 the 
function V\(ki,k2,ks) by VA(k F (ki),k F (k 2 ),k F (k 3 )) = 
V\(ki,k 2 ,ks) for all wave vectors fcj in the same patch 
ki, where i = 1,2, 3. 

VI. PARAMETERS 

The initial condition for the flow of the couplings is 
given by Hubbard interactions Va (k\, k 2 , fc 3 ) = U. For 
most results discussed here we take U — 3t. We choose 
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this rather strong initial interaction because we are inter- 
ested in the breakdown of the Landau-Fermi liquid due 
to interaction effects and do not aim at a classification of 
possible weak coupling instabilities. For the most results 
shown here t' — 0.3t, which is in the range of the values 
reported for the cuprates. 

We will vary mainly the temperature T and the par- 
ticle density near and below half filling via the chemical 
potential fx. We considered values between \i = — 0.7t 
and fi = — 1.35t which corresponds to fillings between 
fa 99% and « 62% of half-filling. The van Hove filling, 
where the FS exactly touches the saddle points, is given 
by fi = -At' = -1.2t. 

For a given \i the dependence of the average particle 
number on T is weak and irrelevant for the results. Both 
parameters \i and T change the effective phase space for 
the various scattering processes. In particular, increased 
temperature provides a larger phase space for particle- 
hole processes with momentum transfer (7r,7r), which 
play an important role. This is similar to the quasi- ID 
organic conductors where above a certain temperature 
the band curvature due to interchain hopping becomes 
irrelevant and—ID nesting effects determine the behavior 
of the systemE3. 

Typically we start at energy scale Ao = it and inte- 
grate the flow down to the critical scale A c which we 
define as the scale where the first coupling reaches a high 
value like 50i. 



VII. COUPLINGS AND SUSCEPTIBILITIES 

We will discuss the results of our numerical RG scheme 
by analyzing the interaction on the FS and susceptibili- 
ties. In the flow to strong coupling we will identify the 
most relevant, i.e. divergent couplings. In the absence 
of scale invariance we cannot expect to obtain simple 
expressions for the form of their divergence, therefore we 
use their numerical values as function of the scale to make 
a qualitative comparison. 

Along with the interactions we calculate the d-wave 
pairing susceptibility Xdw for zero pair momentum and 
the spin susceptibility Xa(o) around = (tt, 7r). The 
method we use is described elsewherecj. Typically both 
Xdw and Xs(<t) grow strongly as the interactions flow to 
strong coupling. In general the ratio of these suscepti- 
bilities is a complicated, non-monotonic, function of the 
scale. Therefore we do not attempt to draw sharp bound- 
aries between different cases. 

As discussed earlier in Sec. I it is often useful to analyze 
also the response the coupling to uniform external charge 
and spin fields given by 

Hc/S = J J^y hc ^ ( 4,t c m ± c L c *u)' 

where the subscripts c and s stand for charge and spin 
respectively. The fc-independent bare couplings h e , a de- 



velop a fc-dependence due to vertex renormalizations de- 
pending on k leading to dressed vertex functions h c / s {k). 
We cannot directly incorporate these renormalizations of 
the uniform external fields in our present RG scheme 
with IR cut-off as they involve only excitations in a low 
energy region of width T around the FS. Therefore we 
calculate h c / s (k) using the RPA for the effective theory 
below the IR cut-off A. This means that we use the scale- 
dependent interactions V\(ki, k%, k 3 ) in the one- loop ver- 
tex corrections for the external couplings, which then can 
be summed up and solved for h c / s (k,A) (see Appendix 
A). 

VIII. RESULTS: THREE REGIMES 

In the density range we examined, we always found a 
flow towards strong coupling at sufficiently low temper- 
ature. The character of this divergence of the coupling 
functions varies continuously with density and temper- 
ature. However we can identify three qualitatively dif- 
ferent regimes, which we will call the d-wave dominated 
regime, the saddle point regime and the approximate nest- 
ing regime, as illustrated in Fig. Our analysis does 
not allow us to draw sharp boundaries between the differ- 
ent regions. Rather, the character of the strong coupling 
flow changes in a crossover-like fashion as one moves from 
one region into the other. 

In order to show the main features we examine the flow 
for three densities typical for each regime. The Fermi 
surfaces and locations of the patch centers that label our 
coupling constants are displayed in Fig. ||. In Figs. 
and H we show snapshots of the couplings at the scale 
where the largest coupling have exceeded the order of 
the bandwidth: we plot the dependence of the coupling 
V\(ki, k2,k 3 ) with the first outgoing wave vector fc 3 fixed 
at point 1 closest to the saddle points or at point 3 closer 
to the BZ diagonal. In Fig. ^ we compare the flow of 
several relevant couplings as a function the RG scale and 
in Fig. [l^ we plot the behavior of the d-wave pairing 
susceptibility Xdw and the AF susceptibility x s (7r, 7r). In 
the following we describe the three regimes in detail. 
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FIG. 6. Fermi surfaces and the 32 points for the three dif- 
ferent chemical potentials discussed in the text, (n) denotes 
the average particle number per site, i.e. (n) = 1 corresponds 
to half-filling. The dots on the FS (solid line) indicate the 
patch centers with patch indices given by the numbers. The 
dashed line denotes the Umklapp surface (US). 
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a. The d-wave dominated regime: At band fillings 
around the van Hove filling fi = —1.2t and low tempera- 
ture T = O.Olt (see right plots in Figs. 0, |, and |), the 
divergence of the coupling functions only occurs at a low 
scale and the d-wave pair scatterings are by far the most 
strongly divergent couplings. In Fig. [?], they appear as 
red and blue features along the lines with patch num- 
bers \ki — k%\ — N/2 (on these lines, the incoming pair 
momentum is zero). Other couplings like Umklapp and 
forward scatterings (the red, violet and black lines in Fig. 
^) grow, too, but are much smaller than the Cooper cou- 
plings. This is the typical flow to strong coupling in the 
lightly shaded regions in Fig. At low temperatures 
it extends also to densities slightly higher than the van 
Hove density. 

A closer analysis shows that the d-wave component 
in the pair scattering is generated at intermediate scales 
by the particle-hole processes with momentum transfer 
(tt, 7r) corresponding to the second term in Eq. (|). This 
type of flow to strong coupling can be considered as a 
Kohn-Luttinger-type Cooper instability where the repul- 
sive scattering in the particle-hole between the saddle 
points first generates a sizable initial value for the d- 
wave pair scattering and is then gradually cut off at lower 
scales because the phase space for the (it, 7r)-particle-hole 
processes decreases due to the shape of the FS. 

The dominance of the d-wave Cooper scattering is also 
seen in the comparison of the susceptibilities: the <i-wave 
pairing susceptibility Xdw grows much faster than the AF 
susceptibility Xs{o) ( see Fig. |To| ) . 

The uniform charge susceptibility k is somewhat sup- 
pressed at intermediate scales but very close to the in- 
stability the attractive Cooper scatterings in the forward 
scattering channel start to dominate the vertex correc- 
tions to the charge coupling and cause a pole in the RPA- 



like expression (see Eq. Bl ) for h c (k) for k near the saddle 
points. This is then the reason of a sharp upturn in n 
(see Fig. [ll]) as also observed by Halboth et al£3. At low 
scales the uniform spin susceptibility Xs(0) is suppressed 
to zero by the strong attractive g 4 couplings favoring sin- 
glet formation. However at higher scales, which are not 
related to the flow to strong coupling, the naive Stoner 
criterion for ferromagnetism is fulfilled due the large DOS 
around the van Hove filling. As discussed in the intro- 
duction, we ignore this effect. 

b. The saddle point regime: Next we increase the 
temperature to T = 0.04i and choose a band filling 
slightly above the van Hove filling such that the FS 
crosses the US (chemical potential fi — —t). Now the 
scale where the couplings reach the order of the band- 
width is strongly enhanced. In Fig. we observe that 
next to the d-wave pair scatterings new features have 
developed. The strongly repulsive interactions, for in- 
stance (fci,fc 2 ) ~ (24,25) — > (fc 3 ,fc 4 ) w (1,17), corre- 
spond to (73-type Umklapp scatterings which now diverge 
together with the repulsive Cooper couplings. The for- 



ward scatterings of 52-type also show a strong increase 
towards the divergence. In addition there is a general in- 
crease for couplings with momentum transfer (jv, it) due 
to the enhanced influence of the particle-hole channel 
with this momentum transfer. On the other hand we 
also observe strongly attractive couplings emerging e.g. 
(fci,fc 2 ) ~ (16,17) — ► (k 3 ,k 4 ) w (1,18). These processes 
correspond to Umklapp 94 processes of pairs with both 
incoming particles at the same saddle-point and outgo- 
ing particles on opposite sides of the FS. Since these 
pairs have small total momentum they couple into the 
Cooper channel and are driven to strong attraction along 
with the attractive Cooper couplings with zero pair mo- 
mentum. This clearly demonstrates that Umklapp and 
Cooper channel are strongly coupled. For this choice of 
parameters the AF susceptibility grows considerably to- 
wards the divergence and is as large as the d-wave pairing 
susceptibility. 

We call this the saddle point regime because the flow 
to strong coupling is dominated by the saddle point re- 
gions. Here, as we will show, we find the key signatures 
of the ISL and the basic mechanism of the two-patch 
model described in Sect. Ill is at work: the diverging g%- 
type Umklapp scattering between the saddle point re- 
gions drives the forward scattering of 52-type to strong 
repulsion, correspondingly the coupling h c {k) of external 
charge fields to these FS parts and thus their contribu- 
tion to the charge compressibility k is increasingly sup- 
pressed as we approach the instability. This can be seen 
and r 



in Figs. 
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L2[ In contrast to the FS near the saddle 
points regions, which tends to incompressible i.e. insulat- 
ing behavior, for k in the BZ diagonal the charge coupling 
h c {k) is more or less unchanged. Therefore one can ex- 
pect that at wave vectors near (w/2,tt/2) gapless charge 
excitations will remain while near the saddle points the 
FS will be truncated. We use this particular behavior of 
the fc-space local charge compressibility to define the sad- 
dle point regime (darker gray regions in Fig. 13): Here 
the charge couplings around the saddle points continue 
to go to zero if we integrate the flow far out of the pertur- 
bative range without any indication of the upturn in k 
that we found in the ci-wave dominated regime. However 
we repeat that in this analysis the border between sad- 
dle point and d-wave dominated regime is a continuous 
cross-over. 

The uniform spin susceptibility exhibits a similar al- 
beit somewhat more isotropic suppression (see Figs. [H] 
and |l2|) when we approach the instability in this sad- 
dle point regime. On the one hand this is plausible be- 
cause the rapid growth of the d-wave susceptibility sig- 
nals strong singlet pairing tendencies. On the other hand 
the AF susceptibility XsC 71 ") 71 ") seems to diverge as well, 
from which one might expect long range AF order, i.e. a 
strong coupling state with nonzero Xs(0)- 
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FIG. 7. (color) Snapshot of the couplings Va(&i, ^3) with first outgoing wave vector £3 fixed at point 1 (see Fig. ^) 



when the largest couplings have exceeded the order of the bandwidth for the three different choices of chemical potential and 
temperature discussed in the text. The colorbars indicate the values of the couplings. 
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FIG. 8. (color) Snapshot of the couplings Va(/ci, hi, fe) with first outgoing wave vector hz fixed at point 3 (see Fig. ^) 
when the largest couplings have exceeded the order of the bandwidth for the three different choices of chemical potential and 
temperature discussed in the text. The colorbars indicate the values of the couplings. For hi = 22, hi — hi « (7r,7r) for hi, ks. 
close to the US. 
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FIG. 9. (color) Flow of the couplings for 32-patch system: d-wave Cooper (blue dashed lines), 33 Umklapp (e.g. 
14(24,24,1),^ (23, 23, 2), solid red lines), g 2 forward (e.g. Va(24, 1, 24), Va(23, 2, 23), black dashed dotted lines), g 4 Umk- 
lapp couplings (e.g. V^(16, 17, 1), solid violet lines) and Umklapp scatterings Va(21,21,4) in the BZ diagonal (green) for the 
three different choices of chemical potential and temperature discussed in the text. 
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Here we argue that for the saddle point regime the 
more likely candidate is a spin liquid state with strong 
short range AF correlation but nonzero spin gap. As 
explained, the flow to strong coupling in this regime is 
caused by the coupling and mutual reinforcement of the 
rf-wave pairing and the Umklapp processes between the 
broad saddle point regions. Therefore the strong cou- 
pling state should feature the singlet pairing of the d-wave 
channel and a strong enhancement of \s (o) for q w (jr, it) . 
This is exactly what we observe for Xdw & n d Xs (o) ■ More- 
over due to the extension of the saddle point regions the 
peak of Xs(<t) is very broad around (tt,tt) and does not 
sharpen significantly in the flow. Therefore we expect 
a rather short AF correlation length of 2-3 lattice spac- 
ings. This is in contrast to the t' = case very close to 
half filling, where we find sharp peaks in Xs(o) developing 
around q — (7r, 7r) and where one would expect AF long 
range order at T = 0. 

c. The approximate nesting regime: The plots on the 
left in Figs. ^, |[ and || show the flow for a higher filling 
(/j = — 0.8t). In this case the leading interactions are 
Umklapp couplings between the BZ regions where the 
FS intersects the US (see red features in Fig. ||) while 
the importance of the vicinity of the saddle points de- 
creases. We call this the approximate nesting regime. 
The dominating FS regions are now further away from 
the saddle points due to the higher filling. As a conse- 
quence the coupling between Umklapp and pairing chan- 
nel decreases and the d-wave pairing processes become 
less relevant. This can be seen best from the weaker flow 
of the attractive Cooper couplings in Fig. ||. Now the 
AF susceptibility clearly exceeds the d-wave pairing sus- 
ceptibility. This signals increasing AF ordering tenden- 
cies which are in accordance with sharper (-7T, 7r) features 
in the interactions (see Fig. decreasing suppression 
of Xs(0) relative to its initial value (see Fig. 11), and a 
sharper peak of Xs(o) around (7r,7r). The charge suscep- 
tibility is also suppressed like in the saddle point regime, 
however the FS regions with smallest charge couplings 
stay fixed to the US and therefore move towards the BZ 
diagonal if we increase the filling. 

We emphasize that in our RG treatment the next near- 
est neighbor hopping t' is important for the existence of 
a sizable saddle point regime. For zero or very small t' 
the FS is closer to the US in the BZ diagonals and the 
(it, 7r) scattering between the rather flat FS faces domi- 
nates even more strongly than in our approximate nesting 
regime with more FS curvature. If we now decrease the 
bandpfiJling, at some point, as pointed out by Zheleznyak 
et al.ES aneLexplicitly shown for the 2D case by Zanchi 
and SchulzO, these processes are cut off at low scales and 
can only serve as generators of an attractive d-wave initial 
condition. With t' very small the system crosses rather 
sharply from a nesting regime into a d-wave dominated 
regime without going through a saddle point regime in 
between. 



H=-1.2t, T=0.01t 




0.5 "0 0.5 "0 0.5 

A/t A/t A/t 

FIG. 10. d-wave (heavy solid line) and AF susceptibil- 
ity(heavy dashed line) for the three different choices of chemi- 
cal potential and temperature discussed in the text. The thin 
lines denote the flow of the bare susceptibilities without ver- 
tex corrections. The mark at the A-axis indicates the scale, 
where the largest coupling reaches lOt. 
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FIG. 11. a) Flow of the charge compressibility k normal- 
ized to their value at the initial scale Aq = 4t for n = —1.2t 
(dashed dotted line), fi = —t (solid line) and n = — 0.8t 
(dashed line), b) Flow of the uniform spin susceptibilities 
normalized to their initial values for fj, = —t (solid line), 
fi = — 0.8t (dashed line), [i = — 0.6t (dashed-dotted line), 
and [i = — OAt (dotted line). For increasing electron density, 
Xs(0, A)/xs(0, Ao) is less suppressed at low scales. For all 
curves T = 0.04i. 




FIG. 12. Change of the charge h c (k) (left plot a)) and spin 
couplings h s (k) (b)) normalized to their initial values of quasi- 
particles with wave vector k on the FS as the electronic inter- 
actions flow to strong coupling from a 96 point calculation at 
fi = —t and T — 0.04t. The different lines are for points close 
to the saddle points (solid lines), and points closer to the BZ 
diagonal (dashed lines). 
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FIG. 13. Dependence of the flow to strong coupling on the 
chemical potential fi and temperature T for t' = 0.3i and ini- 
tial interaction U = 3t. Above the thick line we can integrate 
the flow down to zero scale without reaching an instability. 
Below the thin broken line the d-wave pairing susceptibility 
Xdw exceeds the AF susceptibility Xs( 7r , 71 ") when the largest 
couplings have reached the order of the band width. Above 
this line, Xs(K,n) is larger than Xdw- The darker gray region 
denotes the saddle point regime where the charge coupling 
of the saddle point regions goes to zero and the total charge 
compressibility is suppressed. The lightly shaded region rep- 
resents the d-wave dominated regime. Left to the thick verti- 
cal line the instability is increasingly dominated by couplings 
away from the saddle points, we refer to this region as the 
approximate nesting regime. 

Apart from the suppression of the total charge compress- 
ibility described above there are other potential instabil- 
ities in the forward scattering channel For example, as 
pointed out by Halboth and MetzncrlI3, there appears to 
be a strong tendency towards Labbe-Friedel or Pomer- 
anchuck FS deformations which break the square sym- 
metry. These are mainly rectangular deformation modes 
which split the degeneracy of the saddle-points. However 
we will ignore those tendencies for the reasons discussed 
above. We have checked that a moderate deformation 
of the FS that breaks the square symmetry and leads to 
saddle point splittings of the order of the critical scale 
~ O.lt does not invalidate the results described above. 

A difference to the two-patch analysis of Sec. Ill is that 
the saddle point regime, where we observe the ISL sig- 
natures in our iV-patch calculation, is found at positive 
temperatures and densities slightly higher than the van 
Hove density assumed in the two-patch analysis. The rea- 
son for the latter is that in the A^-patch flow the FS parts 
away from the saddle points reinforce mainly the Cooper 
channel. Only if the FS really crosses the US there is suf- 
ficient low energy phase space for the Umklapp processes 
which then act together with the Cooper processes lead- 
ing to the unusual strong coupling flow. For similar rea- 
sons nonzero temperature is needed for the saddle point 



regime. Moderate T smears out the FS and provides 
additional phase space for both the particle-particle pro- 
cesses with small total momentum and the particle-hole 
processes with momentum transfer (jr, 7r). Especially due 
to the latter there is a certain temperature range where 
this thermal phase space gain outweighs the ordinary de- 
crease of the one-loop contributions for increasing T and 
the critical scale A c is enhanced with respect to its T = 
value. 



IX. DISCUSSION AND CONCLUSIONS 

We have presented an A^-patch renormalization group 
analysis of the 2D Hubbard model and found indications 
that the path from a Fermi liquid-like state to the Mott 
insulating state may pass through a spin liquid phase 
with partially truncated FS and incompressible regions 
around the saddle points. Certainly the above results 
have to be interpreted with care and are only qualitative 
as they are an attempt to learn about possible strong 
coupling states from extrapolating weak-coupling flows. 
However they demonstrate that the breakdown of a Fermi 
liquid through an ISL state with partially truncated FS 
seems to be a viable concept, because in the saddle point 
regime the qualitative features of the ISL, e.g. spin and 
charge gaps, are visible as tendencies in our weak cou- 
pling approach. The essential phenomenon which can be 
identified as the cause for the ISL in the two-patch model, 
namely the coupling of Umklapp and pairing channel, is 
found to exist in a sizable temperature and density range 
also in our improved RG calculation, which includes the 
entire Fermi surface. We believe that this behaviour is 
robust because it only requires sufficiently large low en- 
ergy phase space around the saddle points but does not 
rely on further details of the interaction or dispersion 
relation. What remains to be clarified is when this in- 
terplay between pairing and Umklapp processes, which 
frustrates symmetry-breaking tendencies and thus leads 
to an ISL, indeed represents an energetically favorable 
situation for the system. Another interesting and related 
aspect is the question of the precise conditions for which 
the overlap between the channels becomes too small, such 
that at T — the system can still undergo a transition 
into a symmetry-broken state with presumably renormal- 
ized properties. In our calculation such symmetry-broken 
states are suggested on either side of the saddle point 
regime, e.g. in the d-wave dominated phase or closer to 
half-filling in the approximate nesting regime. 

Our approach certainly bears some appealing features 
when compared to the high-T c cuprates. However note 
that especially very close to half filling, in the approx- 
imate nesting regime, our description will be much too 
simple, as interaction effects which are not taken into ac- 
count will become large. On the other hand further away 
from half filling in the saddle point regime we can hope 
to give a reasonable qualitative description of the driv- 
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ing forces for the breakdown of the Landau-Fermi liquid. 
Due to the mutual reinforcement between Cooper and 
Umklapp channel d-wave pairing correlations appear in 
a natural way at an enhanced scale on the threshold to 
the Mott state. If the insulating tendencies are strong 
enough they will lead to the ISL formation around the 
saddle points. 

The stabilization of the ISL in the vicinity of the sad- 
dle points opens up a novel channel to enhance Cooper 
pairing on the remaining open parts of the FS. A simi- 
laij-*nechanism was recently proposed by Geshkenbein et 
al.Ej who examined a model with infinite mass preformed 
pairs existing at higher temperatures also in the vicinity 
of the saddle points. Let us assume that the ISL has 
formed in a region (called the A region) around the sad- 
dle points at an energy scale Aisl- Then the dominant 
coupling between the ISL and the open FS parts (called 
the B regions) will occur through the exchange of zero 
momentum hole pairs in the Cooper channel. Further it 
will occur in the d-wave pairing channel. We denote by 
Vab the pair scattering matrix element between the ISL 
in the ^-regions and the open B regions of the FS at the 
scale Aisl and by ca the energy relative to the chemical 
potential to add a hole pair to the ISL. Then, at energy 
scales A < Aisl, an additional attraction Vbb is gener- 
ated between pairs in the open B regions, which has a 
pole at 



A c = Aisl exp 



Nai 



n A 



Nab 
dk 



(9) 



Vl B {k) 



B-FS 



(2tt) 2 VF {k) 



Here the integral is over the Fermi surface of the B re- 
gions and tia denotes the number of intermediate states 
with two additional particles in the A regions per lattice 
site. 

Although we do not have a full theory of the strong 
coupling phase, it is plausible to assume that in the A- 
regions a charge gap spreads out along the US in anal- 
ogy to ladder systems which when lightly doped show [Si- 
multaneously channels with and without charge gaptrEM 
Then the open S-rcgion of the FS will enclose an area 
measured from the US, determined by the hole density 
and so the superfluid density will be given by the hole 
density in the saddle point regime. Also there will be 
two energy scales, a higher one determining the onset of 
the ISL in the A regions, and a lower scale setting the 
transition temperature to superconductivity, T c . These 
features are in nice qualitative agreement with the ob- 
servations in the underdoped cuprates. However a full 
microscopic theory of the strong coupling phase and also 
the crossover to the more conventional d-wave supercon- 
ductivity in the lower electron density d-wave regime re- 
mains to be worked out. Note that in the latter regime 
Umklapp scattering is irrelevant at low energy scales and 
the superfluid density is determined by the electron den- 
sity, not the hole density. 



Finally we note that the ISL concept might provide a 
microscopic basis for understanding the ARPES results 
on the cuprates which clearly-show the truncation of the 
FS around the-saddle pointsEJ and also for phenomeno- 
logical modelsEa which have proven to be plausible de- 
scriptions of the transport properties of the normal state 
of the underdoped high-T c cuprates. 



ACKNOWLEDGMENTS 

We would like to thank G. Blatter, J. Feldman, D. 
Geshkenbein, L. Ioffe, H. Knorrer, W. Kohn, U. Leder- 
mann, K. LeHur, W. Metzner, M. Troyer, E. Trubowitz, 
F.C. Zhang, and M. Zhitomirsky for valuable discussions. 
C.H. acknowledges financial support by the Swiss Na- 
tional Fonds. 



APPENDIX A: THE ONE-LOOP SELF-ENERGY 
AND THE FERMI SURFACE SHIFT 

Here is a short overview of the results for the RG flow of 
the Fermi surface with fixed particle number. In order 
to obtain the FS flow we calculate in every RG step the 
change of the one-loop selfenergy given by the contribu- 
tions in Fig.|[ Due to the approximations made for the 
couplings this selfenergy is constant over a single patch 
and only yields a patch-dependent shift of the Fermi sur- 
face. In order to keep the particle number fixed we adjust 
the chemical potential after each step. Quite generally we 
find that the FS parts which have the strongest repulsive 
scatterings with momentum transfer close to (tt, tt) de- 
velop positive self-energies and are therefore are shifted 
inwards during the RG flow. The reason for that be- 
comes clear if one considers a model interaction which is 
sharply peaked and repulsive at Q = (tt,tt). For the self 
energy E(fc) of a particle with wave vector k on the FS 
one primarily has to examine the Hartree term (first term 
in Fig. [|), which is the main contribution for the typi- 
cal divergence of the couplings. This diagram contains 
a propagator with a differentiated cut-off function and 
gives a positive contribution if the state k + Q is occu- 
pied and zero contribution otherwise. After subtraction 
of the FS average of S(A;f) (or more precisely a constant 
which keeps the particle number fixed) this yields a pos- 
itive selfenergy £(fc) > for particles outside the US 
(because then in general the state k + Q is occupied for 
a FS with the densities and t' values we are interested 
in) and a negative shift E(fe) < for states inside the 
US. In our case the interaction only has a broad peak 
around therefore in general also FS points inside 

the BZ can be pushed inwards provided they are more 
affected by this repulsion than the average FS (this hap- 
pens in the overdoped /i = — 1.3t case). The flow of the 
the selfenergies with fixed FS close to the instability is 
shown in Fig. b) for different positions on the FS: for 
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the FS points near the saddle points T,\(k) flows to pos- 
itive values, while for /cf in the BZ diagonal it becomes 
negative. The resulting movement of the FS points if 
we include T,\(k) in the dispersion, i.e. allow the FS to 
move, can be seen in Fig.|lJ a). It reveals the tendency 
of the FS to become flat- thus remaining in the vicinity 
of the Umklapp surfaceEa. In both cases the density is 
kept fixed at (n) w 0.88 per site. Our RG results are 
in qualitative agreement with calculations using a model 
interaction due-to AF spin fluctuations!^!] and the FLEX 
appr oximat ior£j . 




FIG. 14. Left: Initial (open circles) and final (squares) FS 
for t' = O.St and (n) as 0.88 per site. Right: Flow of the 
self-energy on the FS (solid line: point closest to the saddle 
points, dotted line: point closest to BZ diagonal). The flow 
was stopped when the largest coupling reached the bandwidth 
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APPENDIX B: CALCULATION OF THE 
UNIFORM SUSCEPTIBILITIES 



|l5| . The kernel $(fc) is the to = 0, and then q — > limit 
of the Matsubara sum of the product of two propagators, 
and given by the derivative of the Fermi function: 

G>(k\ — lim n F(t(k+q))-riF(e(k)) _ dnp\ _ / R r,\ 

1 ' ~ g™0 e(k+q)-e(k)) ~ dE \ E =e(k)- \P*) 



k+q 

v<-<--<K>-<C 

FIG. 15. Diagrammatic expression for the renormalization 
of the couplings h c / s (k) to external charge or magnetic fields. 
For uniform external fields we take q — > 0. 

The uniform susceptibilities are then given as 

«=-/ -^h° c (k)Hk)h c (k), (B3) 
Xs(0) = "/ ^ h °s(kWk)h s (k). (B4) 

In absence of an instability the coupling functions for 
zero momentum transfer V\(k, k', k) — \V A {k, k', k') and 
— kV\{k, k' , k') would converge to the Landau interaction 

functions f s (k, k') and f a (k, k'), respectively, and the ex- 
pressions for the susceptibilities obtained with the above 
scheme reduce to the results from Fermi liquid theory. 



The uniform (q — > 0) susceptibilities describing the re- 
sponse to external charges and magnetic fields cannot 
be calculated successively by lowering the IR cut-off as 
they only involve degrees of freedom very close to the FS 
(the width of this region is given by the temperature). 
Therefore we determine these responses for the effective 
theory below the cut-off A with the interactions at this 
scale as the effective interactions renormalizing the cou- 
pling to the external fields via vertex corrections. More 
precisely we calculate the effective couplings hi(k) [i = c 
for charge and i = s for spin) of quasiparticles on the FS, 
occurring in the Hamiltonian as 

- — r7rh r /„(k) (ct cr * ± ct cr , ). 
(27r) 2 ' M fe >T k,l k 'l' 

Denoting the bare coupling as h?^ s we can express the 
effective coupling as 

-> ( dk 1 -» -> -> -> 
hi(k) = h°{k) + I -j—- hi(k)$(k')Vi{k, k') i = c,s 



(Bl) 



where V c (k,k') = (-V A (k,k',k') + 2V A (k,k',k)) for the 
charge and V s (k, k') = — V\(k, k', k') for the spin cou- 
pling. Diagrammatically this equation is shown in Fig. 
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